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Within the previously developed Dubna-Mainz-Taipei meson-exchange model, the singularity 
structure of the tvN scattering amplitudes has been investigated. For all partial waves up to F waves 
and cm. energies up to Ty ~ 2 GeV, the T-matrix poles have been calculated by three different 
^^ ■ techniques: analytic continuation into the complex energy plane, speed-plot and the regularization 

CN ' method. For all 4-star resonances, we find a perfect agreement between the analytic continuation 

and the regularization method. We also find resonance poles for resonances that are not so well 
^ , established, but in these cases the pole positions and residues obtained by analytic continua- 

tion can substantially differ from the results predicted by the speed-plot and regularization methods. 
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"q . I. INTRODUCTION 



Ever since the A(1232) resonance was discovered by Fermi and collaborators in 1952 [l|-l3|, the excitation spectrum 
of the nucleon has played a fundamental role in our understanding of low-energy hadronic physics. The most direct 
*Zl . evidence for resonance structure is based on pion-nucleon elastic and charge-exchange scattering. Because total 
0\| , angular momentum, parity, and isospin are conserved within the realm of the strong interaction, the S matrix for 
_i , the reactions tt + N — > n' + N' may be decomposed into the partial wave amplitudes T/j_, with / the isospin, 
r — I . i the orbital angular momentum, and the ± indicating the total spin of the hadronic system, J = i zt ■^. Two 
^^ , decades after the discovery of the A(1232), a dedicated program at the meson factories had provided enough data to 
establish a rich resonance spectrum of the nucleon [J, |5| . The further partial wave analysis was driven by studies 



cd 



of the Karlsruhe-Helsinki (KH) [y, lZ[ and Carnegie-Mellon-Berkeley (CMB) [8|, |9[ collaborations. In the following 
k> , years, R. Arndt and coworkers at VPI compiled the data base SAID [10|, which was later updated and extended in 
^ ! collaboration with GWU [llj, |l2[. The works of the CMB, KH, and GW/VPI groups are the main sources of the 
nucleon resonance listings in the Review of Particle Physics (PDG) |13| . 

In the most intuitive way, a resonance is an intermediate state of target and projectile that lives longer than in a 
typical scattering process. Translated into the language of scattering theory, resonances are defined as poles of the 
S matrix on unphysical Riemann sheets. Different methods were developed to derive the resonance properties from 
the observables. In the 1930's it was suggested that a Breit-Wigner function should be a good representation for a 
resonance pole, and the Breit-Wigner formula for spin zero particles and its generalization to finite spin were devel- 
oped (see an illustrative discussion in Cottingham and Greenwood, p. 241 in Ref. llJ]). Later the discussion centered 
on the rapid increase of the eigenphase shifts through 90° and on the related backward looping of Argand diagrams [l5| . 

However, Breit-Wigner parametrizations have been found to be very model dependent. As was recently shown in the 



framework of effective quantum field theory, Breit-Wigner masses are in general field-redefinition dependent [16| . The 
same model dependency also applies to electromagnetic properties as charges, magnetic moments, transition moments, 
and form factors. On the other side, all these resonance properties are uniquely defined at the pole of the 5" matrix [17| . 

The analytic properties of the S matrix are imposed by the principles of unitarity and causality. Because of 
unitarity, each physical channel leads to a square-root branch point of the partial wave amplitude T{W) at the 
respective threshold, with the result that T{W) is a multi-valued function in the complex W plane. In particular, a 
partial wave for elastic nN scattering is described by the amplitudes T^^l on the physical and T^^l on the unphysical 
sheet. The experimental amplitudes are identified with the amplitudes above the cut, Texp(VF) = lime_j.o T^^l (W + ie): 
with W > M + m and e > 0. As a consequence, the physical sheet has a discontinuity over the real axis, along the 
right-hand cut M -\- m < W < oo, with M the nucleon and m the pion mass. In addition, a relativistic theory has 
a left-hand cut starting at M — m. Causality requires that the physical sheet be free of any further singularity, the 
nucleon resonances should appear as simple poles on the unphysical sheet closest to the real axis of the physical sheet, 
in agreement with Hohler's remark 18|: "It is 'noncontroversial among theorists' (see Chew \1 w and the references 



in my 'pole-emics', p. 697 in Ref. \2(\ 1) that in S-matrix theory the effects of resonances follow from first order poles 



in the 2nd sheet." A pole on the second sheet, described by T'-'^'{W) « rp/{Mp ~W — iTp/2), will often lead to a 
maximum of the experimental cross section near W = Mp. The resonance is therefore defined by (i) its pole position 
in the complex cm. energy plane at Wp ~ Mp — i Tp/2, with Mp the real part of the pole position and Tp the width 
of the resonance, and (ii) the residue of the amplitude, r^ ~ \rp\cxp{i 9p) , at the pole. 

The focus of the present work is on how to extract resonance properties from the data, that is, how do we find a 
pole in the complex energy plane having at our disposal only data on the real axis and use different pole-extraction 
methods to extract one, and only one pole per resonance. The problem is how to eliminate the background from 
the experimental amplitude and to isolate the pole contribution. The time delay method and the related speed plot 
technique turned out to be the most quantitative tools to detect poles and extract their parameters from the data. 
A resonance is characterized by a maximum time delay, the time passing between the arrival of a wave packet and 
its departure from the collision region. In general, a large time delay indicates the formation of an unstable particle 
in the intermediate state. However, misleading effects can occur by rapid variation of backgrounds, such as narrow 
cusp effects above 5'-wave thresholds or spurious singularities due to phenomenological parametrization of form 
factors and cut-offs. If a resonance lives long enough, it should decay into all energetically possible final states, unless 
prevented by general selection rules. Furthermore, the pole position derived from the data should not depend on how 
the resonance is excited or decays. Whereas such resonances exist in atomic and also in nuclear physics, some caveat 
is in order for nucleon resonances. As an example, a simple classical model of the A(1232) leads to the conclusion 
that the pion stays in its orbit around the nucleon for only about 100° of a full circle. 
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The pole positions of nucleon resonances have been traditionally derived from the data by the speed plot 
which is related to the time delay. The "speed" is defined by the slope of the amplitude with energy, dT{W)/dW, 
which eliminates constant backgrounds. The speed plot shows the modulus of the speed, SP(W^) — \dT{W)/dW\, 
and resonances are identified with peaks in the plot. The resonance parameters are then obtained by fitting the speed 
of a single pole to the data at physical values of the energy. The idea of the speed plot has been recently generalized 



to higher derivatives by the "regularization method" (RM) [22|. Within a convergence circle given by the closest 
neighboring singularity, the Laurent expansion about a pole is given by the sum of T^°^°{W) and a Taylor series 



T'^^^{W). With an increasing number of differentiations, the signature of the pole sticks out more and more sharply, 
whereas more and more leading terms of the Taylor series disappear. It goes without saying that the differentiation 
of the data will fail after a few steps because of numerical instabilities. However, the regularization method is an 
interesting tool to study the singularity structure of analytic models. In particular, this method will reveal any rapid 
variation due to cusp effects and spurious singularities introduced by phenomenological parameterizations. 



The inherent model dependence of all the methods has been often overlooked. Examples are the decomposition 
of the scattering amplitude in a Breit-Wigner form and a background [20| and the fact that eigenphases of a 
multi-channel system usually do not run through 90° because of the Neumann- Wigner no-crossing theorem [23|. 
With regard to the latter point, Dalitz and Moorhouse jl5| stated: "With such a complexity of branch cuts without 
physical significance, we must conclude that the eigenphase representation for the S matrix is not generally a useful 
representation for the scattering in the neighborhood of a resonance" . On the other hand, a single-channel procedure 
can not reveal much information about resonances 'far away' from that channel, that is, if the particular resonance 
has only a small branching ratio for the particular channel. 



The dual nature of the nucleon resonances, the "bare" resonance based on quark configurations and the "dressed" 
pion-nucleon continuum state, has been studied by several collaborations within Lagrangian models, in particular 
by EBAC at Jefferson Lab 2J, |25[, the groups at Giessen 26[, Jiilich 27 1, Valencia 28 1, and Darmstadt |29| as 
well as the Dubna-Mainz-Taipei collaboration [30j . Such approaches are also the basis for structure investigations of 
resonances by photo- and electroproduction. 



Although Lattice QCD (LQCD) has obtained promising results for the masses and several resonances of the 



low-lying hadrons J3ll433| . the large pion mass used and the quenching ap prox imation make it (yet) impossible to 



treat the resonances as a pion-nucleon scattering state. As shown by Ref. |3J] for the A^A form factors, the chiral 
extrapolation from the stable A at large pion masses to the experimental pion mass yields unexpected and rather 
dramatic non-analytic effects at the A — > Ntt threshold. The very fact that lattice theory can not yet describe the 
pion-nucleon final state interaction, makes it impossible to compare the lattice data to the experimental scattering 
amplitudes in a direct way. However, the LQCD data could be compared to the results of a dynamical model obtained 
with the same pion mass as used in LQCD. As long as the pion mass is large, m,^ > tojv* — rn^, the N* appears 
as an excited bound state. In this case the mass predicted by LQCD should be compared to the "dressed" mass of 
the dynamical model, which contains the "bare" mass and a real self-energy due to meson loops. In the mass region 
where the N* decays, the LQCD calculation would provide a complex amplitude, comparable to the scattering phases 
of a partial wave expansion. Based on the work of Liischer [35|, the width of the rho meson [36[ and of the Delta 
resonance |37| have been recently studied. In such a case the lattice data can be treated like experimental amplitudes. 



that is, by speed plot or Breit-Wigner analyses. However, as long as the LQCD pion mass lies above the physical value, 
an extrapolation to the physical pion mass will still be necessary by use of a dynamical model or an effective field theory. 



In the present contribution we study the Dubna-Mainz-Taipei model (DMT) 30| by comparing the pole parameters 
resulting from analytic continuation with approximate procedures such as the speed plot and the regularization 
method. The DMT model is a field-theoretical meson-exchange model for ttN scattering. In contrast to the partial 
wave amplitudes obtained in Refs. [8l-ll2l|. which inevitably involve a smoothing of the data, the energy dependence of 
the DMT amplitudes is largely determined from theoretical considerations even though there are free parameters in 



the model. Another strong point of the DMT model is the fact that the prescribed background contribution, a major 
source of the model dependence for the resonance extraction procedure, has been found to give excellent agreement 



if applied to low energy electromagnetic pion production [38|- In Sect. II we give an overview of the speed- plot, 
time-delay, and regularization methods to derive the pole parameters. The Dubna-Mainz- Taipei model is presented 
in Sect. Ill, in particular with regard to the definitions of resonant vs. background terms as well as form factors and 
cut-offs. Our results for the resonance parameters are reported in Sect. IV, and the different techniques to derive 
these parameters are compared. We conclude with a summary and outlook in Sect. V. 



II. HOW TO FIND A RESONANCE 
A. Time delay 

In the framework of a potential scattering problem, resonance phenomena are related to the formation and the decay 
of intermediate quasi-stationary states. An ideal resonance should decay in all energetically possible final states, unless 
forbidden by selection rules for a specific channel. Provided that the interaction is of short range, resonances can 
be characterized by the time delay between the arrival of a wave packet and its departure from the collision region. 
In general, the time delay shows a pronounced peak at the resonance energies. The time delay was introduced by 
Eisenbud in his Ph.D. thesis [39| and later applied to multi-channel scattering theory by Smith [40|. Following the 
work of Refs. 



15 



41 



43| . we define the time delay for a single-channel scattering problem in a partial wave as follows: 
, , , / 1 dS(W)\ dS(W) 

where S{W) = exp[2iS{W)] is the S matrix and 6{W) the scattering phase shift. A simple ansatz for a unitary S 
matrix with a pole at Wp = Mp — ^Tp and a constant background phase Sb is given by 

S(W) = ^^P + i^P - ^ ^2rSs ^ ^2^S^iW) 2rSu (2) 

where 5ii{W) = arctan[irp/(Mp — W)] is the resonant phase. The related T-matrix T = [S — l)/(2i) and the real 
matrix K = i{\ — S){1 + S)^^ take the forms 

n»')-]S;r|^«»'-+sWe»«., (3) 

^ ' Mp-W - \Tptan5B 

Combining Eqs. ([1]) and ([2|), we obtain a simple Brcit-Wigner form for the time delay. 

In this ideal case, the maximum time delay is At{Mp) = 4/rp. We observe that both the real (Mp) and imaginary 
(— iPp) parts of the pole position are determined by the time delay (and therefore also by the scattering matrix) at 
real (physical!) values of the cm. energy W. 



In order to describe a system of coupled channels, the lifetime matrix Q was introduced 



H, 



,m--^'-^suw). (6) 



The analog of the time delay for a multi-channel system was found to be the trace of Q as function of W . This trace 



takes a more transparent form after diagonalization of the S matrix 15|, |44 1 



"[owi=2i:^, (7) 

a 

with 5a the eigenphase shifts. However, a realistic application of the lifetime matrix requires the knowledge of all open 
channels, that is, the reactions irN — > nN , ttN — > ttttN, ttttN — > ttttN, ttN — > r/N, and so on. As a consequence of 
Wigner's no-crossing theorem [45|, individual eigenphases have a complicated energy dependence. It is therefore only 
the sum of the eigenphases that shows distinct resonance structures. Based on this observation, it has been recently 
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proposed to search for resonance parameters by studying the traces of multi-channel T and K matrices [46|, 
respective scattering matrices were constructed from experimental data for the ttN and rjN channels and models for 
the two-pion channels. 

B. Speed plot 



The speed plot of a partial wave amplitude T is defined by 

dT{W) 



SP{W) 



dW 



(8) 



As was recognized by the Particle Data Group already in the early 1970's, the speed plot is a convenient tool to 
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extract the pole position of a resonance [5|, |48|. This technique was intensively studied by Hohler who wanted to 



extract resonance parameters from the Karlsruhe-Helsinki partial wave analysis (KH80) [2l|, |49|. Because the KH 
analysis was restricted to elastic pion-nucleon scattering, a multi-channel treatment like the construction of the 
lifetime matrix was out of reach. 

For a single channel system, the time delay and the speed plot are identical up to a factor 2. In particular, the 
ansatz of Eq. ([3]) leads to the speed 

As a result, the speed plot shows a maximum for W = Mp, which defines the real part of the pole position in the 
complex W plane. The imaginary part of the pole position can be obtained from the relation 

SP(A/p±irp) = iSP(Mp). (10) 

In practical applications, this straightforward method is numerically not very stable and fails completely for the 
5ii(1535) resonance, which is very close to the t] threshold. In order to allow for inelasticities, we therefore assume 
the functional form 



nw) = „, J ,^ (11) 

with Tp, the complex residue at the pole, given by 



Mp ~W -\Tp 



^^^T{W)\y,^^^^-\Tp\e'^. (12) 



For a single (elastic) channel, the comparison with Eq. ([3]) yields the residue. The more general ansatz Eq. (|TT|) leads 
to the speed 

which is fitted to speed data obtained from the partial wave amplitudes in the vicinity of the maximum. The pole 
parameters Mp^Tp and \rp\ are then obtained from the best fit to the selected speed points. Finally, the phase 0p of 
the residue is obtained from the phase of dT/dW, 



Im (dT/dW) 
tanei^^ ^ ' ^ 



^ Re (dT/dW) 



(14) 



C. Regularization method 

As discussed before, the speed plot technique is successful if the background under the resonance is approximately 
constant. It fails if the background changes over the resonance region. The regularization method extends the idea 



behind the speed plot by construction of higher derivatives, T*^^) = d'^ T/dT^ , N > 1 2^- Let there be an analytic 
function T{z) of a complex variable z with a first-order pole at some complex point /.j = .t + iy. This function can be 
any of the T-matrix elements, and the variable z is identified with the cm. energy W in order to compare with the 
speed plot technique. The described function takes the form 



T{z) = + T{z) , (15) 




smooth background 



where fi and r are the position and residue of the pole. Of course, the experiment can determine the T-matrix 
elements only for real values of W. In order to continue T(W) into the complex energy plane and to search for the 
pole position, we construct a regular function / by multiplying T with the factor ij, — z, 

f{z) = {fi-z)T{z), (16) 

with f{fJ.) = r. In the neighborhood of the pole, the function / can be expanded in a Taylor series. Because the 
scattering matrix can be accessed for real arguments only, we construct /(/i) from the derivatives of / taken along 
the real axis, 

71 = 

This expansion is explicitly written to order A^, and Rn{W, fj,) stands for the higher orders. The derivatives /'"■'(VF) 
can be turned in derivatives of the T matrix by use of Eq. (|16p , and mathematical induction leads to the following 
equation: 

/(")(VK) = {^i-W) r(")(VF) - n T^"-^\W). (18) 

Insertion of these derivatives into Eq. ([16]) cancels all the terms in the sum except for the last one, 

" = ^^""^J-^^ (m - M/)(^+^) + Rn{W, /.) . (19) 



In the neighborhood of the pole, the remainder R^ should decrease with increasing N. Assuming that the higher 
derivatives can be neglected for a sufficiently large value of N and taking the absolute values of both sides, we obtain 
an approximation of the pole parameters at 0{N), 

On condition that the Taylor series converges and in the limit A^ — >■ oo, rjsr and /iw should approach the values r 
and /x, respectively. In the next step, we (i) write the pole position as a general complex number, /x = a + ib, (ii) 
raise both sides of Eq. ([^0]) to the power of 2/{N + 1), and collect the information about the T-matrix and the pole 
position on the right and left side, respectively. The result is a parabolic equation in W, 



{gn - W)^ + 9 



N_ _ jv + i 



-^^ • (21) 



This equation relates the pole position (a = Mp, b = — ^ Tp) and the absolute value of the residue, \r\, to the T-matrix 
values on the real axis, as obtained from a model or an energy-dependent partial-wave analysis of the data. Finally, 
the phase of the residue is determined by 

ImT(^)(VF) 



tanWjv = 



RcT(^){W) 



(22) 

W=M„ 



The comparison of the above equations with the results of Sect. IIIBI shows that the speed plot is identical to the 
regularization method for A^ = 1. 

The further procedure is as follows: (i) Construct the A'*'' derivative of the T-matrix element and the right-hand 
side of Eq. ([2T|) . Note that the pole parameters are uniquely determined by the exact knowledge of T{W) in only 
three points. However, the problem is to choose the right points. If the distance between the points gets too large, 
the influence of other singularities may increase. If the points are too close, numerical problems may occur, (ii) 
Solve Eq. (|2T|) for the pole parameters by either choosing various three-point sets to evaluate the right-hand side and 
perform a statistic analysis of the results or fitting the right-hand side of the equation to a three-parameter parabolic 
function. In our approach we have chosen the latter option. 

In closing this section we note that the regularization method does not depend on any particular functional form of 
the T matrix. However, we have to assume that (a) the A'*'' derivative can be constructed with a sufficient precision 
and (b) the pole position lies within the circle of convergence for the Taylor expansion, that is, no further singularities 
should intrude into the region between the pole and the related resonance region on the real W axis. 

D. Poles from analytic continuation 

The most accurate way to determine pole positions and residues is certainly obtained by analytic continuation into 
the complex region. Because resonance poles can not appear on the physical sheet, we have to take a careful look at 
the structure of different Riemann sheets opening at all branch points for particle production in a coupled-channels 
model. The most important particle thresholds in our energy region below 2 GeV are ttttN, ttA, tjN, and pN with 
branch points at 1178 MeV, (1350 - 50z) MeV, 1486 MeV, and (1713 - 75i) MeV, in order. In the dynamical DMT 
model we have included the ttN and mrN channels in all the partial waves and the r]N channel in the Su partial 



wave. However, the -kttN channel is treated in a phenomenological way as will be described in the following section. 
Because the particular ansatz for the two-pion width, Eq. (|43|). contains only even powers of q2TT, the model does not 
give additional branch points for the two-pion channel. This leads to a relatively simple sheet structure, and we can 
easily reach the poles on the most important second Riemann sheet (first unphysical sheet). Technically, we first map 
the relevant region on this sheet by contour plots and search for the accurate pole position by standard root finding 
routines applied to the function 

M^.) = ^ = 0. (23) 

Next we obtain the residue by approaching the pole along different paths in the complex plane, 

Resnz)l^=\m,{z-Zp)T{z). (24) 

A word of caution is in order. If the fitted form factor parameters, e.g., Aq of Eq. ([57|) or Xji of Eq. (P5|) become 
smaller than about 500 MeV, additional poles can appear in the region where the resonance poles are expected. In 
order to avoid such spurious singularities, it is very important to map out the structure of the T matrix very precisely. 
Also the speed-plot and regularization methods are very helpful to distinguish between resonance and spurious poles, 
because the latter ones usually affect T matrix at real VF in a similar way as a broad background does. 



III. THE DUBNA-MAINZ- TAIPEI MESON-EXCHANGE MODEL 

The Dubna-Mainz-Ta ipe i (D MT) ttN meson-exchange model was developed on the basis of the Taipei- Argonne ttN 



meson-exchange model 50|-|52| which describes pion-nucleon scattering up to 400 MeV pion laboratory ener gy. The 
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DMT model was extended to cm. energies W = 2.Q GeV by inclusion of higher resonances and the r]N channel 

The model gives an excellent fit to both nN phase shifts and inelasticity parameters in all the channels up to the 

F waves and energies of 2 GeV, except for the Fn partial wave. The DMT -kN model is also a main ingredient of 



the Dubna-Mainz-Taipei dynamical model describing the photo- and clectroproduction of pions |54| up to 2 GeV. In 
particular, thi s model gives an excellent agreement with pion production data from threshold to the first resonance 
region 



38 



55| . In this section, we briefiy outline the ingredients of the Taipei- Argonne ttTV model and then describe 



some relevant features of the DMT meson-exchange coupled-channels model. 

A. Taipei-Argonne meson-exchange iiN model 

The Taipei-Argonne model describes the elastic scattering of pions and nuclcons. It is based on a three-dimensional 
reduction of the Bethe-Salpeter equation for an effective Lagrangian involving tt, N ^ A, p, and a fields. Let us first 
describe some kinematics for -kN scattering, 

^{q)+N[p)^Ti{q') + N{p'), (25) 

where q,p,q', and p' are the four- momenta of the respective particles. The total and relative four- momenta are 
P = p + q and k ~ pil-nis) — (7?77v(s), respectively, where s ~ P^ ~ W"^ is the Mandelstam variable. The dimensionless 
variables »7^(s) and r]N{s) represent the freedom in choosing a three-dimensional reduction, they are constrained by 
the condition rjN + '77r = 1 • For further details we refer to Ref. 



The Bethe-Salpeter equation for ttN scattering takes the form 

TttN — SttTV + Bt^nGoTttN , (26) 

where B^^n is the sum of all irreducible two-particle Feynman amplitudes and Gq the free relativistic pion-nucleon 
propagator. Equation (1^51) can be cast into the form 

Tttn = Bttn + B-^nGoT-^n , (27) 

with 

Bt:N = -BttAT + Bt,n{Go — Go)i?7rJV , (28) 

where Go(fc; P) is an appropriate propagator to obtain a three-dimensional reduction of Eq. (|26p . This propagator is 
chosen to maintain the two-body unitarity by reproducing the ttN clastic cut. However, there is still a wide range of 



possible propagators satisfying this constraint. A standard choice is the Cooper- Jennings reduction scheme 56 1 



1 f ds' 

(2^ 

X 6^+\[vn{s')P' + k]^ - m%) <5(+)([7^(s')P' - k]^ - ml) , (29) 



Goik;P) = j—-^ / —J{s,s')[a{s,s')r+^ + mN 



with P' = \J s' / s P and the superscript (+) of the (5- functions standing for the positive energy part of the propagator. 
The variables / and a are dimensionless variables containing the freedom of reduction, they are constrained by the 
conditions /(s, s) = 1 and a(s, s) ~ ?7jv(s) to ensure the reproduction of the elastic cut. Furthermore, Cooper and 
Jennings define [56| 



a{s,s)=riN{s), f{s,s)^— -^ ^^2 ' (^O) 

ss' — [mf^ — rn^Y 

which leads to the following expression for Go in the cm. frame: 

(2^)3 V^-ys- Vi+ys-^^'^^ 4i?jv(fc)K(fc) 

where EN{k) and ETj{k) are the nucleon and pion energies for the three-momentum k, J^ ~ _Ejv(fc) + ET,{k) is the 
total energy in the cm. frame, and ri{s, k) = E]^{k) — rip^{s^),/s^ . With these relations we obtain the following ttN 
scattering equation: 



t{k', k; W) = v{k', k;W)+ / dk"v{k', k"; W)gQ{k"; W)t{k" , k; W) . (32) 
The explicit relations between the quantities defined by Eqs. (|26l) and ([32]) are given by 

t{k',k;W) = f dk'„dknS{k'„-fj')T{k',k;W)S{kn-fi), (33) 

v{k', k;W) ^ I dk'„dkoS{k'o - f/)B{k', k; W)S{ko - f/), (34) 

goik;W) = fdkoGoik;W), (35) 

with fj' = fi{sg, k') and f/ = fjisj:, k). 



10 
The effective Lagrangian involving the Tr,N,a,p, and A(1232) fields takes the form 

^(0) _ (v) 

+ ^A^[g^'^ - {Z +^)n'']fANN ■ d^Tf , (36) 

with A^ the Rarita-Schwinger field operator for the A resonance and Taw the isospin transition operator between 
nucleon and A. The driving term BttN in Eq. p6p is approximated by the tree diagrams of the interaction La- 
grangian of Eq. p6p , the direct and crossed N and A diagrams as well as the t-channel a- and p-exchange contributions. 



Furthermore, the driving term v{k, k'; W) of Eq. (j34p is regularized by covariant form factors associated with each 
leg of the vertices, 



Kipl) 



(37) 



with Pa the four-momentum and rua the mass of particle a. In the present work we choose the value n = 10 as used 
by Ref. [53]. 



As in the work of Afnan and collaborators |57l.l5a|. the Pn phase shift is constrained by imposing the nucleon pole 
condition. This treatment leads to a proper renormalization of both nucleon mass and ttNN coupling constant. It 
also yields the necessary cancelation between the repulsive nucleon pole contribution and the attractive background, 
such that a reasonable fit to the TriV phase shifts in the Pn channel can be achieved. 

The following parameters were allowed to vary in the fit procedure: the products gaNNga-w-n, gaNNga-Jn, and 
gpNNgpiTiT as well as 5 for the t-channel a and p exchanges, tti^ ' 9ttNA' ^^'^ ^ ^'^^ ^^^ ^ mechanism, and the cut-off 
parameters Aq, of the form factors given by Eq. ([37]). The experimental ttN phase shifts were well describe d up 



to pion laboratory energies of 400 MeV. The resulting parameters and predicted phase shifts can be found in Ref. [52 1 



B. DMT meson-exchange model including 77 A'^ channel and higher resonances 

As the energy gets larger, heavier mass channels like aN, rjN, ttA, and pN as well as a non-resonant continuum 
of t:t:N states become increasingly important, and at the same time more and more nucleon resonances appear as 
intermediate states. For each contributing resonance R, the Hilbcrt space is enlarged by a bare resonance R which 
acquires a width by its coupling to the ttN and r]N channels through the Lagrangian 

^x - igi^RRrN ■ TT + igl^i^^RNv + h.c. , (38) 

where N, R, tt, and 77 denote the field operators for the nucleon, bare resonance i?, pion, and eta meson, respectively. 
The full t-matrix can be written as a system of coupled equations, 

Uj [W) = V,, {W) + Y, v,k [W] gk (W) tk, iW) , (39) 

k 
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with i and j denoting the tt and r] channels and W is the total cm. energy. 

In general, the potential Vij is the sum of non-resonant (vf^) and bare resonance (vf^j) terms, 

v,j{W)^vf^{W)+v^^{W). (40) 

The non-resonant term v^^ for the ttN elastic channel, as defined in Sec. IIII Al contains contributions from s- and 
M-ehanncl Born terms as well as i-ehannel contributions with to, p, and a exchange. The parameters in v^^-are fixed 
by an analysis of the pion scattering phase shifts for the S and P waves at energies W < 1300 MeV 5^1 ■ In the 
channels involving the rj, the potential u^^ is assumed to vanish because of the small 7]NN coupling [59|. 

The bare resonance contribution arises from the excitation and decay of the resonance R, 

-5W = ^^^^' (41) 

where Mj^ and hi J denote the bare mass for the resonance R and the bare vertex operator R -^ irf-q + N, respectively. 
The matrix elements of the potential v^{W) is symbolically expressed by 



MA,,q;W)gi°^gf^f,{A,,q';W) 



v^iM. ?'; W) = ^'^'^"^"^^^(„)^^7:^'^"•^ , (42) 



W - M^^^ + ^T^j^ [W) 



where q and q' are the pion (or eta) momenta in the initial and final states, and g^j. denotes the resonance vertex 
couplings. As in Ref. [S^l, we associate a covariant form factor with each leg of the vertices. Therefore, each vertex 
fi in Eq. P^ contains three form factors as given by Eq. ((57)) . and Aj stands for a triple of cut-offs, (Ajv, A/{, A^). In 
Eq. (|^ we have also included a phenomenological term T'^j^[W) in the resonance propagator to account for the -kttN 
decay channel. Therefore, our "bare" resonance propagator already contains a phenomenological "dressing" effect due 
to the coupling to the 'ktiN channel. With this prescription we assume that any further non-resonant coupling to the 
ttttN channel can be neglected. Following Refs. 60|, |61| we parameterize the two-pion width by 

ru xi+qi ^'-' 

qo J \XJ, + qh 



rlTW^rrM^) 4f^] ' (43) 



where i is the pion orbital momentum and (727r = q27riW) the momentum of the compound two-pion system. Further- 
more, the two-pion width and two-pion momentum at resonance are denoted by F^ and qo, respectively. We note 
that this form accounts for the correct energy behavior of the phase space near the three-body threshold [60| . In our 
present work, Tj^ and X^ are considered as free parameters. Therefore, each resonance is generally described by 
6 free parameters, the bare mass M^', the decay width F^'^^ ', two bare coupling constants gl ' and 5^ , and two 
cut-off parameters Aji and Xj^. The generalization of the coupled-channel model to the case of N resonances with 
the same quantum numbers is then given by 

N 

vf^{q,q';W) = Y,vff{q,q';W), (44) 

n=l 

with 6 free parameters for each resonance. 

The solutions of the coupled-channel equations of Eq. (p9)) . with potentials given in Eqs. (|40ll43p . were fitted to the 
experimental irN phase shifts and inelasticities by variation of the bare resonance parameters. The fit gave a good 
agreement with the data for all channels up to the F waves and energies below 2 GeV, except for the partial wave 

Fn 



3 of . The predictions of the DMT model for the resonance parameters are presented in the next section. 
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IV. RESULTS AND DISCUSSION 

In this section we present the nucleon resonance parameters as derived from the DMT model. The hsted resonances 
fulfill the following criteria: (i) the pole position is restricted by Mp < 2 GeV and Tp < 0.4 GcV, (ii) the residue 
is larger than about 1 MeV, and (iii) the branching ratio for the one-pion channel is limited by 2rp/rp > 10 %. 
Furthermore, the pole position obtained by the renormalization method has to be stable over a range of derivatives 
(N values). 

A. S- waves 



As reported in Ref. [53[, we need four Sn resonances to fit the ttN scattering amplitude in this channel, instead of 
the three resonances listed by the PDG [13|. The additional resonance 5'ii(1878) was found to play an important role 
in pion photoproduction as well |53|. but was not seen in both the ttN — > r/N reaction and recent measurements of rj 



photoproduction from the proton [62l l63j. However, in the analysis of Durand et al. j64l |65| on eta photoproduction 
a new Su state of mass M = 1707 MeV and width F = 222 MeV was needed in order to get good agreement with 
the data. 

In our analysis (see Fig. [H left panel) we find four poles, three of them below 2 GeV. The exact pole positions are 
shown by asterisks with a size proportional to their relative strength. Furthermore, the ranges of the PDG pole values 
are displayed by open boxes. The first resonance S'ii(1535)**** is well found by both speed-plot and regularization 
methods. The second state, the 5ii(1650)**** is even better described, with some improvement by the regularization 
method. The third state is very weakly excited in the ttN interaction and can not be seen by the SP, whereas the 
RM finds a state close by. A fourth state, slightly above 2 GeV and classified by PDG as 1-star 511(2090)* is also 
found by SP and RM at close-by positions. 

For the isospin 3/2 state we obtain three poles in agreement with PDG. The first one, 5'3i(1620)**** is nicely 
described by SP and perfectly by RM. The second pole, 5*31(1900)** is similarly well seen by both SP and RGM. The 
third resonance, 53i(2150)*, lies far down in the imaginary region. Only the SP locates a state in the expected energy 
region, however, with only a small fraction of the residue can be obtained. 

B. P-viraves 

For the P-wave resonances we show our results in Fig. [2] In the Pn partial wave we get 3 poles, all of them being 
somehow related to the 3 PDG states. However, our second and third pole positions lie further down in the negative 
imaginary region, the third pole is even below the range shown in the figure. The position of the first resonance, the 
Roper Pii(1440)**** is very well described by both SP and RM. The two higher states have very weak signals in the 
ttN amplitude. Although the RM yields a considerable improvement over the SP, it also misses the exact pole position. 

From the results in Table |ll] we find that only a very small fraction of the residue can be found with these 
techniques. This is a general property for poles lying far away from the real axis. A detailed analysis of the contour 
for this partial wave also shows a zero position sitting closely above the pole position, and therefore masking the pole 
if searched for along the real axis. 



13 



In the Pi3 partial wave we also find three poles. Two of them are very well reproduced by RM. PDG lists only two 
states, the Pi3(1720)**** with a large error bar for the imaginary part and the Pi3(1900)**, however, with no pole 
position data given. Our results for the Pi3(1720)**** lie close to the PDG values, with the imaginary part close to 
the lower limit of PDG error bar. 

The isospin-3/2 partial wave, P31 is an exceptional case. Here the lowest state is not the dominant one, and also 
the PDG lists a P3i(1750)** without pole values. We find a lowest pole near 1800 MeV with a very large imaginary 
value outside of our considered range. The dominant P3i(1900)**** lies only slightly outside the PDG box and is 
very well described by RM. 

The Delta resonance, ^33(1282)****, has of course the largest strength of all the poles found in our analysis. It is 
very well described by both the SP and RM techniques and coincides with PDG. Two further, rather weakly excited 
P33 states arc found at higher energies, which can be related to the states P33(1600)*** and P33(1920)*** reported by 
PDG. While the analytic result for the P33(1600)*** agrees fairly well with RM, the P33(1920)*** lies outside of the 
considered range. 

C. D- waves 

In Fig. |3]we show the poles of D-wave resonances. PDG lists three J = 3/2 states with isospin 1/2, Z?i3(1520)****, 
-Di3(1700)*** and Di3(2080)**, but even a total number of four regions with extracted pole values, three of these 
PDG boxes show up in our figure. Our analysis agrees very well within the reported ranges, where we also locate 
four resonance states. While the dominant first resonance is perfectly described by SP and RM, the SP can only 
locate another (third) pole. The improved RM technique, however, is able to find all four positions with high accuracy. 

For the D15 partial wave we find a rather simple contour, and 2 pole positions can be well located in perfect 
agreement with SP and RM. However, only the 1)15(1675)**** lies in our considered range. 

For the isospin-3/2 D-wave resonances, PDG reports only one 4-star resonance, the Z?33(1700)***. Our analytic 
results are reproduced well by SP and perfectly by RM. A similar agreement is also found for the resonances 1)33(1940)* 
and 1)35(11930)***. All 3 resonances are located very close to the PDG boxes. 

D. F- waves 

Figure |4] displays the P-wavc resonances. The Pi5(1680)*** is the most important resonance in the 3'''^ resonance 
region and perfectly described by SP and RM. We also find a weak second F15 state. However, because of its 
location close to the real axis, this state is also well described by SP and RM. In the vicinity, PDG reports the 
Pi5(2000)** resonance, however, with no pole position given. In the Fn partial wave we can not find a resonance 
by our analysis, whereas PDG lists a Pi7(1990)**. Our result is based on the fact that the inclusion of a bare reso- 
nance in this wave does not significantly improve the x^ fit to the data. A similar conclusion was found in Refs. |lll.ll2|. 

For the isospin-3/2 resonances, we find three poles as also reported by PDG. The first poles in both partial waves. 
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i^35(1905)*** and i^37(1950)***, arc well described by SP and RM. However, our positions are significantly above the 
PDG bounds situated at large imaginary values. This disagreement is not too surprising, because the uncertainty of 
a pole position rises with the distance from the real axis. We also find a third pole close to the ^35(2000)**, which is 
however very weak and not well seen in the physical region. 

V. SUMMARY AND CONCLUSION 

Within our previously developed Dubna-Mainz-Taipei (DMT) dynamical model we have studied the pole structure 
of the pion-nucleon T matrix for all S*, P, D and F partial waves in the energy range up to cm. energy W = 2.0 GeV. 
We have analytically continued the T matrix to the unphysical region in the second Riemann sheet and have 
accurately calculated the pole positions and the residues down to IvaW > —200 MeV. These exact values were 
compared with the results obtained by the speed plot, based on the first derivative of the T matrix, and a recently 
developed regularization method using higher derivatives of T. 

Our general conclusions are: 

• The number and positions of the DMT poles, obtained by the newly developed analytic continuation method, 
are in very good agreement with the current PDG results [13[. 

• All 4-star resonances are nicely reproduced by both the speed-plot technique and the renormalization method. In 
many cases, however, the latter approach gives significant improvements for the resonance parameters, especially 
for the pole residues. 

• We have confirmed the need for a forth Sn resonance in addition to the three resonances listed by PDG [1 
The additional 5*11(1878) resonance was found to play an important role in pion photoproduction 



not seen in both the iiN — > rjN reaction and recent data of rj photoproduction from the proton [62 



531, but was 



• We have shown that the regularization method is a reliable method to extract the pole parameters from single- 
channel data. In the absence of a full experimental knowledge about all the channels, this method can be 
sequentially applied in order to track down the poles relevant for the experimentally known channels of a 
multichannel amplitude. 
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FIG. 1: The pole positions for the S waves in the complex energy plane. The (red) stars show the results found by analytic 
continuation, the (green) solid circles and the (blue) open circles are determined by the speed-plot and the regularization meth- 
ods, respectively. The rectangular regions show the range of the pole positions listed by the Particle Data Group (PDG08) [1^ ). 
The size of the (red) stars is proportional to jr-pl/Fp, and therefore a measure for the strength of the resonance poles. 
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FIG. 2: The pole positions for the P waves in the complex energy plane. The notation is the same as in Fig. [T] 
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FIG. 3: The pole positions for the D waves in the complex energy plane. The notation is the same as in Fig. [T] 
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TABLE I: Pole positions Wp = Mp 



^iTp and absolute values of the residues 



at the pole, all in MeV, as well as the 



phases 9p of the residues for S-wave resonances. The first lines give the exact pole positions and residues as obtained by 
analytic continuation of the DMT model, the second and third lines show the values obtained by the speed-plot (SP) and the 
regularization method (RM(N)) with N being the largest value for stable derivatives. The star classification and the numerical 
values of the PDG are listed in the last line. 



N* 


ReWp 


-ImWp 


1 rp 1 


Opideg) 


511(1535) 


1499 


39 


14 


-45 


SP 


1499 


26 


7 


-47 


RM(1) 


1499 


26 


7 


-47 


+*** 


1510 ± 20 


85 ±40 


96 ±63 


15 ±45 


5ii(1650) 


1631 


60 


35 


-83 


SP 


1642 


49 


22 


-74 


RM(6) 


1631 


52 


28 


-119 


**** 


1655 ± 15 


83 ±8 


55 ±15 


-75 ±25 


6-11(1880) 


1733 


90 


16 


-29 


SP 


- 


- 


- 


- 


RM(6) 


1806 


83 


10 


-172 


new 










531 (1620) 


1598 


74 


23 


-98 


SP 


1598 


69 


23 


-99 


RM(6) 


1598 


75 


24 


-100 


**** 


1600 ± 10 


59 ±2 


16 ±3 


-110 ±20 


531 (1900) 


1774 


36 


3.8 


179 


SP 


1775 


18 


1.0 


-166 


RM(5) 


1777 


28 


1 


-157 


** 


1870 ± 40 


90 ±25 


10 ±3 


-20 ± 40 
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TABLE II: Pole positions residues for P-wave resonances. For further notation see Tab. [T] 



N* 


ReWp 


—ImWp 


kp| 


9p{deg) 


Pii(1440) 


1371 


95 


50 


-79 


SP 


1366 


90 


48 


-87 


RM(5) 


1371 


95 


50 


-78 


**** 


1365 ± 15 


95 ±15 


46 ±10 


-100 ±35 


Pii(1710) 


1746 


184 


11 


-54 


SP 


1721 


92 


5 


-164 


RM(6) 


1756 


150 


11 


-49 


*** 


1720 ± 50 


115 ±75 


10 ±4 


-175 ±35 


Pi3(1720) 


1693 


136 


20 


-43 


SP 


1683 


118 


15 


-64 


RM(4) 


1695 


133 


19 


-34 


**** 


1675 ± 15 


98 ±40 


13 ±7 


-139 ±51 


P3i(1910) 


1900 


87 


13 


-116 


SP 


1896 


68 


7 


-118 


RM(6) 


1901 


87 


10 


-113 


+ !|!** 


1880 ± 30 


100 ± 20 


20 ±4 


-90 ± 30 


^33(1232) 


1212 


49 


49 


-42 


SP 


1218 


44 


41 


-35 


RM(2) 


1218 


45 


41 


-35 


+ !|!** 


1210 ± 1 


50 ±1 


53 ±2 


-47 ±1 


^33(1600) 


1544 


95 


14 


-111 


SP 


1509 


99 


25 


-197 


RM(5) 


1546 


80 


11 


-116 


*** 


1600 ± 100 


150 ± 50 


17 ±4 


-150 ±30 
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TABLE III: Pole positions and residues for D-wave resonances. For further notation see Tab. H] 



N* 


ReWp 


—ImWp 


kp| 


dpideg) 


Di3(1520) 


1515 


60 


40 


-7 


SP 


1516 


61 


40 


-6 


RM(5) 


1516 


60 


40 


-5 


**** 


1510 ±5 


55 ±5 


35 ±3 


-10 ±4 


Oi3(1700) 


1718 


48 


2.8 


-91 


SP 


- 


- 


- 


- 


RM(7) 


1718 


49 


2.9 


-91 


*+* 


1680 ± 50 


50 ±25 


6±3 


0±50 


^13(2080) 


1854 


108 


16 


-97 


SP 


1834 


108 


14 


-134 


RM(6) 


1858 


105 


9 


-83 


** 


1950 ± 170 


100 ± 40 


27 ±22 


~ 


i?i5(1675) 


1657 


66 


24 


-22 


SP 


1657 


66 


24 


-23 


RM(6) 


1657 


67 


25 


-22 


**** 


1660 ± 5 


68 ±6 


29 ±6 


-30 ± 10 


D33(1700) 


1604 


71 


9.4 


-63 


SP 


1609 


67 


9.5 


-52 


RM 


1604 


71 


9.9 


-63 


+*** 


1650 ± 30 


100 ± 20 


13 ±3 


-20 ± 25 


D35(1930) 


1989 


140 


18 


-78 


SP 


1992 


136 


19 


-75 


RM(5) 


1989 


140 


18 


-78 


*** 


1900 ± 50 


133 ± 48 


18 ±6 


-20 ± 40 
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TABLE IV: Pole positions and residues for _F-wave resonances. For further notation see Tab. |I] 



N* 


ReWp 


-ImWp 


kp| 


dpideg) 


Fi5(1680) 


1664 


57 


38 


-26 


SP 


1663 


57 


38 


-28 


RM(6) 


1664 


58 


39 


-26 


**** 


1672 ± 8 


61 ±6 


38 ±3 


-23 ±7 


Fi5(2000) 


1919 


26 


1.0 


15 


SP 


1931 


31 


1.3 


89 


RM(4) 


1919 


27 


1.1 


15 


+* 


1807 


54.5 


60 


-67 


Fi7(1990) 


- 


- 


- 


- 


SP 


- 


- 


- 


- 


RM 


- 


- 


- 


- 


** 


1900 ± 30 


260 ± 60 


9±3 


-60 ± 30 


F35(1905) 


1760 


100 


10 


-66 


SP 


1771 


96 


11 


-47 


RM(5) 


1759 


101 


10 


-66 


**** 


1830 ± 5 


140 ± 10 


25 ±8 


-50 ± 20 


^37(1950) 


1858 


104 


43 


-48 


SP 


1860 


101 


44 


-45 


RM(5) 


1859 


104 


44 


-47 


*+** 


1880 ± 10 


140 ± 10 


50 ±7 


-33 ±8 



